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One of the largest sources of systematic uncertainties in ongoing neutrino-oscillation measurements 
is the description of nuclear effects. Its considerable reduction is expected thanks to the dedicated 
studies of (anti)neutrino-nucleus interactions in the MINERvA experiment. In this article, the 
calculations within the spectral function approach are compared to the charged-current quasielastic 
cross sections reported from MINERvA. The obtained results show that the effect of final-state 
interactions on the (anti)muon kinematics plays a pivotal role in reproducing the experimental 
data. 
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The MINERvA Collaboration has recently reported 
the flux-averaged differential cross section da/dQ^ ec , 
Qi- ec being the reconstructed four-momentum transfer 
squared, for charged-current (CC) quasielastic (QE) scat¬ 
tering of muon antineutrinos 0 on nucleons in the hy¬ 
drocarbon target, CH. In a subsequent paper 0, the 
corresponding result for muon neutrinos has also been 
presented. While five theoretical calculations have been 
compared to the data in Refs. 0,1 , none of them turned 
out to be able to satisfactorily describe both the antineu¬ 
trino and neutrino cross sections. Rather consistent re¬ 
sults have been found using superscaling approaches M 
and the local Fermi gas model with effective interactions 
and the random-phase approximation effects @]. 

In this article, I argue that the effects of final-state 
interactions (FSI) on the ^ kinematics are essential to 
reproduce the MINERvA CC QE data, half of which cor¬ 
respond to Q 2 ec < 0.2 GeV 2 , where the contribution of 
low momentum transfers is significant 0. 

Changing the energy balance in the primary vertex, 
interactions between the struck nucleon and the specta¬ 
tor system affect the final energy of the charged lepton 
E'f.. Additional modifications to E' e , different for i~ and 
£ + , result from an influence of the Coulomb field of the 
nucleus on all charged particles in the final state [g]. 

As known from electron scattering, FSI broaden the 
QE peak and shift it in a kinematics-dependent way, in¬ 
creasing typical Ep at low momentum transfers. To il¬ 
lustrate those effects, in Fig. [Q the predictions of the 
spectral function (SF) approach 0- involving no ad¬ 
justable parameters—are compared to the experimental 
data of Ref. 01. The same mechanism redistributes the 
(anti)neutrino cross sections da/dQ^ ec toward higher val¬ 
ues of Q 2 ec . 

The role of FSI on electron and neutrino cross sections 
in the SF formalism has been studied in detail by the au¬ 
thors of Ref. || for the carbon nucleus. In the approach 
of Ref. IH], the target’s grou nd-state properties are de¬ 
scribed using the hole SF 0 ] , consistently accounting for 


* ankowski@vt.edu 


both the the shell structure 0, 0] and nucleon-nucleon 
correlations [T3 |. To include the effect of Pauli block¬ 
ing, the particle SF calculated following Ref. 0 is em¬ 
ployed. Modification of the struck-particle’s spectrum— 
due to its FSI with the spectator system—is accounted 
for in the energy conservation by making use of the Dirac 
phenomenological potential of Cooper et al. [l5j . FSI- 
induced broadening of the QE peak is described within 
the correlated Glauber approximation 0 , 113 , appl yin g 
the nuclear transparency of carbon reported in Ref. |18[ . 
It is noteworthy that a remarkable agreement with a large 
body of experimental 1 gC(e, e') data 0 . fl9l - l24l ] has been 
observed in Ref. Q. 

Within the approach of Ref. 0], I have calculated the 
CC QE v M and cross sections averaged over the MIN¬ 
ERvA fluxes, expressing both the carbon and free-proton 
(for v lL only) contributions as a function of Q 2 ec defined 
as in Refs. 0 0. I have employed a conservative esti¬ 
mate of the nuclear model’s uncertainties—represented 
by the bands in Fig. [3]— coming from the optical po¬ 
tential parametrization and the contribution of the low- 
excitation energies E x < 26 MeV. Due to the stability 
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FIG. 1. (color online). Double differential C(e, e') cross sec¬ 
tions at scattering angle 13.5° and beam energies 1.35, 1.50, 
and 1.65 GeV. The SF calculations without (dashed line) and 
with (solid lines) FSI effects 3 are compared to the data 





2 



0.01 0.1 1.0 
QL (GeV 2 ) 

FIG. 2. (color online). Differential cross sections da/dQ^ ec 
for CC QE (a) iy, and (b) scattering in MINERvA. The 
SF calculations without (dashed line) and with (solid lines) 
FSI effects [9|] are compared to the data 0 , 0 - The bands 
represent theoretical uncertainties. 


TABLE I. Fit results to the CC QE MINERvA data. 



antineutrino 

neutrino 

combined fit 


including theoretical uncertainties: 

Ma (GeV) 

1.16 ±0.06 

1.17 ±0.06 

1.16 ±0.06 

y 2 /d.o.f. 

0.38 

1.33 

0.93 


neglecting theoretical uncertainties: 

M a (GeV) 

1.15 ±0.10 

1.15 ±0.07 

1.13 ±0.06 

y 2 /d.o.f. 

0.44 

1.38 

1.00 


neglecting FSI (Ma = 

1.16 GeV): 

y 2 /d.o.f. 

2.49 

2.45 

2.42 


of the carbon nucleus, I have assumed that the central 
values of the calculations correspond to E x > 0. To ac¬ 
commodate the contribution of two-body currents in an 
effective manner, the axial mass Ma has been treated as 
a free parameter. 

In calculations of the y 2 distribution, I have accounted 
for correlations between the data, employing the stan¬ 
dard expression 



where the sum runs over all data points, A*, is the differ¬ 
ence between the theoretical and experimental value at 
point k, subject to the total uncertainty 07-, and V de¬ 


notes the correlation matrix reported by the MINERvA 
Collaboration. Including the theoretical uncertainty in 
(jfc, I have added it in quadrature with the experimental 
error. 

The determined Ma values and the corresponding y 2 
per degree of freedom (d.o.f.) are given in Table [TJ The 
difference with respect to Ma = 1-03 GeV, extracted 
predominantly from the deuteron measurements }25l j, in¬ 
dicates a nonvanishing contribution of reaction mecha¬ 
nisms other than those involving a one-body current. In¬ 
terestingly, the effective Ma values turn out to be very 
consistent in the antineutrino and neutrino cases. 

The analyses with and without theoretical uncertain¬ 
ties give rather consistent results for a twofold reason: (i) 
theoretical uncertainties are only relevant at low Q 2 ec , 
where the cross sections exhibit little sensitivity to the 
axial mass value and the calculations reproduce the data 
very well, and (ii) the dominant contribution to y 2 comes 
from the high-(3 2 ec region. 

On the other hand, as shown in TablelH FSI effects turn 
out to be essential to reproduce the experimental cross 
sections of Refs. Hi- It is illustrated also in Fig. [2] 
where the numerical results accounting for FSI and ne¬ 
glecting them are represented by the solid and dashed 
lines, respectively. The axial mass applied in the calcu¬ 
lations is set to 1.16 GeV. 

The employed uncertainty estimate of the nuclear 
model is based on the extensive comparison to electron¬ 
scattering data 0, [H£3 performed in Ref. [§]. It is 
important to note that it conservatively assumes that the 
cross section’s uncertainty is 100% at the kinematics cor¬ 
responding to the excitation energy E x < 26 MeV, where 
the excitations of the giant dipole resonance and bound 
states contribute. Comparisons to low-energy data M- 
|29( |— left for future studies—can be expected to allow this 
assumption to be relaxed, reducing theoretical uncertain¬ 
ties and enabling much clearer discrimination of the cal¬ 
culations with and without FSI than that presented in 
Fig. [2] 

As a final remark, I would like to acknowledge that 
results similar to these presented here have been reported 
from the relativistic Green’s function model [3(|. This 
consistency is likely to be ascribed to the fact that the two 
approaches make use of the same optical potential fl5| . 
Unfortunately, the authors of Ref. [30| have not pointed 
out an essential ingredient bringing their calculations into 
agreement with the MINERvA data. 

In summary, a satisfactory description of the CC QE 
cross sections extracted from MINERvA can be achieved 
when FSI effects on the ^ kinematics are accounted for 
in the nuclear model. It is important to keep in mind that 
the selection of the CC QE signal may exhibit a certain 
dependence on experimental details and is subject to sys¬ 
tematic uncertainties stemming from backgrounds sub¬ 
traction [3l[. Therefore, it is interesting to note that the 
effective values of the axial mass determined in this arti¬ 
cle are in rather good agreement with the result reported 
by the MiniBooNE Collaboration in Ref. j32| . This find- 
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ing can be interpreted as pointing toward a non-negligible 
role of two-body reaction mechanisms at the kinematics 
of the MINERvA experiment. 
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